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LOW-THRUST ORBIT RAISING IN CONTINUOUS SUNLIGHT WHILE THRUSTING 
IN A PLANE PERPENDICULAR TO THE EARTH-SUN LINE 
by James R. Ramler 
Lewis Research Center 

SUMMARY 

An approximate numerical method is derived for determining the motion of an Earth 
satellite following a tightly wound, spiral trajectory which is nearly circular at all times 
while remaining in continuous sunlight. Thrust is provided by electric rockets and is 
continuous, tangential, and in a plane perpendicular to the Earth-Sun line. This differs 
from previous studies which have generally assumed the thrust vector to lie in the orbit 
plane. The maximum time that the satellite can remain in continuous sunlight has been 
computed by optimizing the initial orbit inclination and mission starting time for initial 
thrust-weight ratios from 2. 5x10"® to 8.0x10"® and initial orbit altitudes of 300, 500, and 
800 nautical miles (555. 6 , 926. 0, and 1481. 6 km). Precession of the orbital plane is 
largely accomplished by utilizing the effects of the oblateness of the Earth. The effect of 
reversing the thrust direction during the latter part of the mission on the maxi m um time 
spent in continuous sunlight has been studied. Methods of solution are discussed and re- 
sults are presented as curves. Errors in establishing the required initial orbit conditions 
inherent in any realistic satellite launching are briefly discussed. 


INTRODUCTION 

Testing of electric rockets in an actual space environment will be required to demon- 
strate fully their performance characteristics, in particular, their operation lifetime. 

The operating lifetime of electric rockets may be demonstrated in an actual space en- 
vironment by continuously operating them for a long period of time in orbit about the 
earth. While a number of trajectory concepts may be envisioned, the purpose of this re- 
port is to analyze one of them - that of raising the orbit altitude by thrusting tangentially. 



thus enabling an accurate determination of the thrust level from orbit mechanics. Chang- 
ing the orbit altitude by thrusting tangentially would also demonstrate the capability of 
electric rockets to perform useful missions in the future. 

Such a mission would require continuous electric power to be supplied to the electric 
thrustors for a long period of time. The most likely source of power for this, as well as 
similar missions in the near future, is solar cells which are lightweight and can operate 
for long durations. In order to demonstrate the operating lifetime of electric rockets, it 
will be desirable to keep the satellite with its solar cells in continuous sunlight as long as 
possible. 

For an Earth satellite to remain in continuous sunlight as long as possible, its orbit 
must be near polar and the orbit plane must be rotated with time in an easterly direction. 
This will tend to keep the Sun in a direction normal to the orbit plane, thus continuously 
illuminating the satellite as the Sun changes position in its apparent motion along the 
ecliptic during the year. This orbit rotation or precession can be accomplished by utili- 
zing the effects of the oblateness of the Earth on an inclined orbit (ref. 1). 

The problem of raising the orbit of a satellite using low thrust while remaining in 
continuous sunlight has been studied in reference 2. Altitude raising was accomplished by 
using a continuous tangential thrust in the orbit plane starting from a low circular Earth 
orbit. The maximum final altitude or mission time which could be reached was computed 
by an approximate technique. The effects of initial orbit and mission parameters such as 
inclination, altitude, thrust-weight ratio, and mission starting time were studied. 

Tangential thrusting in the orbit plane has been shown to be an efficient steering pro- 
gram for orbit raising (ref. 3). However, if the solar panels are oriented to always face 
the sun, tangential thrusting in the orbit plane introduces several complications into the 
spacecraft design. For example, the thrustors must be gimbaled with respect to the 
solar panels and a cyclic variation in gimbal angle is required during each orbit revolu- 
tion. The maximum total gimbal angle (equivalent to tj) can become as high as ±70° as 
shown in reference 2 (all symbols are defined in appendix A) . Furthermore, the inertial 
direction of motion of the satellite must be continuously determined onboard the satellite 
in order to thrust in the orbit plane. 

For this study, a steering concept was assumed that would be simpler to imple- 
ment although somewhat less efficient than tangential thrusting. It is based on a concep- 
tual design in which the solar panels and electric thrustors would be rigidly fixed to the 
spacecraft frame (except for possibly a small gimballing capability on the thrustors for 
attitude control) . The solar panel and thrust vector geometry will be discussed in more 
detail later. Briefly, however, the spacecraft would be oriented so that the solar cell 
panels always face the sun. The panels would, therefore, lie in a plane perpendicular 
to the Earth -Sun line. The electric thrustors would lie in the same plane as the solar 
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panels and, in addition, would be oriented to always thrust in the local horizontal plane. 
This geometry will result in a variable component of thrust normal to the orbit plane. 

This study then considered the case of an Earth satellite using electric rockets that 
thrust tangentially in a plane perpendicular to the Earth-Sun line rather than thrusting 
tangentially in the orbit plane as was the case in reference 2 . This steering concept could 
be implemented using only horizon scanners and Sun seekers. 

The thrust vector can be considered to have two components: one in the orbit plane 
and the other normal to the orbit plane. An approximate numerical solution has been de- 
rived to determine the motion of an Earth satellite under these conditions while remaining 
in continuous sunlight. Although this solution was derived to analyze the effect of a nor- 
mal component of thrust resulting from a particular orientation of the satellite with re- 
spect to the orbit plane and the Sun, the approach to the problem and method of solution 
presented in this report would be generally applicable to other low-thrust, circular Earth 
orbit missions having out-of -plane thrusting. The mathematical analysis for this study is 
presented in appendix B. 

By use of the results of appendix B, the maximum time that an Earth satellite can re- 
main in continuous sunlight was computed by optimizing the initial orbit and mission vari- 
ables and by utilizing the oblateness of the Earth as the major orbital-plane rotation me- 
chanism. The ranges of initial variables investigated parametrically included orbit alti- 
tudes of 300, 500, and 800 nautical miles (555.6, 926.0, and 1481.6 km), thrust-weight 
ratios from 2. 5x10” to 8.0x10" and mission starting dates throughout the year. Initial 
orbit inclination and the time of day at the start of the mission were optimized for each 
point. The effect of reversing the thrust of the satellite during the latter part of the mis- 
sion was studied in an attempt to extend the maximum time that the satellite remains in 
continuous sunlight. This reversal would be effected by simply pitching the spacecraft 
over 180°. The time of thrust reversal was optimized for each case. 

For comparison, the case where thrusting is completely tangential and in the orbit 
plane was also investigated, thus extending the results of reference 2 to lower values of 
thrust -weight ratios. 

Errors in establishing the required initial orbit conditions inherent in any realistic 
satellite launching and their effect on the mission are briefly discussed. 


ANALYSIS 

Thrust Vector Geometry 

Consider the spacecraft shown in figure 1 to be in a near polar circular Earth orbit. 
For a satellite orbiting the Earth using continuous low thrust, tangential thrusting has 
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been shown to be nearly the most efficient for increasing energy (ref. 3). As seen in fig- 
ure 1 , tangential thrusting in a circular orbit can be obtained by simply keeping the roll 
axis of the spacecraft parallel to the local horizontal plane. In addition, the solar cell 
panels would be oriented to face the Sun by alining the pitch axis of the spacecraft toward 
the Sun. This will keep the solar cell panels normal to the incident radiation of the Sun, 
thus taking full advantage of their energy conversion efficiency. 

To keep the roll axis of the satellite parallel to the local horizontal plane as it orbits 
the Earth, it would be necessary to control the attitude of the satellite in pitch. Also, to 
keep the pitch axis of the satellite pointed toward the Sun, it would be necessary to con- 
trol the attitude of the satellite in roll and yaw. Rolling of the satellite will have no effect 
on the thrust vector; whereas, yawing will cause it to move out of the orbit plane. With 
the roll axis maintained parallel to the local horizontal plane, however, yawing will only 
result in a component of thrust normal to the orbit plane (the in-plane component remains 
tangential). For this study, then, the thrust vector consists of two variable components; 
one tangential and in the orbit plane and the other normal to the orbit plane. This thrust 
vector alinement may not necessarily be optimum, but, for early research and develop- 
ment flights, it would provide an effective steering program that could be implemented by 
using horizon scanners and Sun seekers. 

With the attitude of the satellite determined as described previously, the thrust 
acceleration vector A will always be tangential and lie parallel to a plane perpendicular 
to the sunline as shown in figure 2. For simplicity, the sunline is defined as a vector 
from the center of the Earth to the center of the Sun, rather than from the spacecraft to 
the Sun; thus, a very small error in angle between the orbit plane and the pitch axis of 
the spacecraft (less than 0. 01° for orbit altitudes up to 10 000 n mi (18 520 km)) is neglec- 
ted. As illustrated in figure 2, attitude changes of the spacecraft can be assumed to be 
periodic over one orbit if the inertial positions of the orbit plane and the Sun are assumed 
not to change during one orbit. This is a good assumption for one or even several orbits 
since the relative motion of the Sun and the orbit plane will be very slow. Because the 
thrustors are fixed parallel to the roll axis of the spacecraft, the thrust acceleration vec- 
tor A will then also periodically change direction with respect to the orbit plane over one 
orbit. As illustrated in figure 2, it can be seen that, in the orbital positions 1 and 3, the 
orbit track and the plane perpendicular to the sunline intersect. When the spacecraft is 
at these positions, the yaw axis will be alined with the radius vector so that there is no 
displacement in roll. However, there will be a yaw angle at positions 1 and 3 equal to the 
angle r] between the orbit plane and the plane perpendicular to the sunline. This will re- 
sult in a component of thrust acceleration W normal to the orbit plane at these positions. 
At positions 2 and 4, which are displaced 90° from positions 1 and 3, there will be no dis- 
placement in yaw but in roll. At these positions, the normal component of acceleration 
will be zero and the total thrust acceleration A will be tangential in the plane of the 
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orbit. From the analysis in appendix B, it is shown that the magnitude of the normal 
component of thrust acceleration W will vary periodically over one orbit from zero at 
positions 2 and 4 to a maximum equal to A sin tj at positions 1 and 3. Likewise, the 
magnitude of the tangential component of thrust acceleration C will vary from a minimum 
equal to A cos 77 at positions 1 and 3 to a maximum equal to the total thrust acceleration 
A at positions 2 and 4. 

The acceleration components, W and C, will also vary secularly with time. Al- 
though it was assumed that the inertial positions of the Sun and the orbit plane do not 
change during one orbit or even during several orbits, they actually do change contin- 
uously with time; that is, the angle rj between the orbit plane and the plane perpendicu- 
lar to the sunline will vary continuously with time. The variation of 77 depends on the 
motion of the Sun in right ascension and declination and the changes in inclination and 
nodal position of the orbit plane. Changes in the inclination and node of the orbit plane 
are caused by perturbations due to the oblateness of the Earth, thrust of the spacecraft, 
and other bodies such as the Sun and Moon. 


Orbit Perturbations 

The forces on a satellite in a circular Earth orbit due to oblateness will cause the 
orbit plane to process about the polar axis of the Earth. An analytical expression for the 
secular rate of change of nodes of an orbit can be obtained by considering only the terms 
through the second harmonic in the potential function of the oblate Earth. The following 
expression is obtained in reference 4 by ignoring the smaller, higher order terms in the 
potential function: 


Jr|(Rj, + cos I for 0°<I<180° (1) 

where J is the coefficient of the second harmonic in the gravitational potential function. 
Note that, for inclinations greater than 90° (retrograde orbits), is positive, which re- 
sults in an eastward precession of the node . This is the desired precession direction to 
keep the Sun in a proper position relative to the orbit plane. All the orbits studied, 
therefore, have inclinations greater than 90°. Perturbations due to the oblateness of the 
Earth of other orbital elements such as inclination, eccentricity, and semimajor axis are 
small and periodic with no resultant secular change. 

The normal component of thrust acceleration W causes secular variations in both 
and I which are given by equations (B12). Other orbital elements will exhibit only 
small periodic changes over one orbit but no secular changes because of the normal com- 
ponent of thrust acceleration. The tangential component of thrust acceleration C changes 
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only the orbit altitude. (Note that, for circular orbits, the tangential- and circumferential- 
thrust vector directions are identical. ) 

As is shown later, perturbations of the orbital elements <fae to the Sun and Moon were 
found to be negligible and are not included in the numerical approximation. 


Earth Shadow Considerations 

By assuming only circular orbits, the geometry of the orbit with respect to illumina- 
tion from the Sun is simplified. The existence of a penumbra region of shadow is neglec- 
ted and it is assumed that the rays of the Sun are parallel, causing a cylinder shaped um- 
bra region behind the Earth. Shown in sketch (a) is an edgewise view of the orbit plane 
and its geometrical relation to the position of the Sun. The angle between the orbit per- 
pendicular (orbital angular -velocity vector) and the sunline is defined as 77. It can be 
seen from sketch (a) that as long as rj is less than 77^, the entire orbit will be in sun- 
light. The situation depicted in sketch (a) is a limiting case in which the orbit is just on 
the edge of the cylindrical shadow and rj equals 77^. 



(The plane of the sketch is the plane containing the sunline and the orbit perpendicular. ) 
From sketch (a), 77^ is defined by 


R 

cos n ^ for 0 ° <Vr>< 90 ° (2) 

«E + h 

This means that, for an orbit to be in continuous sunlight, the orbit perpendicular must 
remain inside a cone having the moving sunline as its axis. The half -angle of this cone 
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is the 7 ]^ corresponding to the instantaneous orbit altitude. This cone (fig. 3, taken from 
ref. 2) is called the cone of tolerance and the intersection of the cone of tolerance and the 
surface of the Earth is called the circle of tolerance. 

An expression for the angle rj can be derived from sketch (b). 



Consider the spherical triangle formed by the orbit plane, the plane perpendicular to the 
sunline, and the equatorial plane. The angle 77 is the angle between the orbit plane and 
the plane perpendicular to the sunline. As seen in figure 3, the angle ^ is defined to be 
the difference between the right ascension of the Sun a and the position of the orbit per- 
pendicular in longitude As such, it is an indication of how far the orbit perpendicu- 

lar is west of (lagging behind) the easterly moving sunline. The orbit index, N = +1, cor- 
responds to the sunline being on the same side of the orbit plane as the orbit perpendicu- 
lar (orbital angular -velocity vector), and N = -1 corresponds to the sunline being on the 
opposite side of the orbit plane as the orbit perpendicular. It should be noted that for 
N = -1, 7 }^ must be 180*^ larger than \f/ for N = +1. 

From sketch (b), for N = -1 

cos(tV + 180°) = 77 + cos I cos(90° - _5)L 

sin I sin(90° - 6) 


Therefore, 


cos 77 = -cos ij/ sin I cos 6 - cos I sin 6 
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For N = +l 


cos xl/ = cos r; + cos I cos(90° + 5) 
sin I sin (90*^ + 6) 


Therefore, 


cos 77 = cos \l/ sin I cos 6 + cos I sin 6 


and, in general. 


cos 77 = N(cos \j/ sin I cos 6 + cos I sin 6 ) (3) 

The angle 77 can then be evaluated from equation (3) and compared with the angle 77 ^ 
evaluated from equation ( 2 ) to determine whether or not the satellite orbit remains in sun- 
light. 


Mechanics of Problem 

During the mission, the orbit perpendicular will describe a path within the circle of 
tolerance. As long as this path remains within the circle of tolerance, the orbit will be in 
continuous sunlight. The problem can then be viewed as starting the mission with the or- 
bit perpendicular inside the circle of tolerance and then maximizing the time that the orbit 
perpendicular remains within the circle of tolerance. Several factors will influence the 
motion of the orbit perpendicular within the circle of tolerance. 

As the Earth moves about the Sun during the year, a and 6 change continuously. 
Although a is an increasing function, 6 varies between approximately 23. 5° and -23. 5°. 
Since the sunline is the axis of the cone of tolerance, as a changes, the cone of toler- 
ance moves in longitude along the ecliptic . As 6 changes, the cone of tolerance moves 
above and below the ecliptic in a periodic manner. Each location of the cone axis corre- 
sponds to a specific time of year. 

Tangential thrust will cause the altitude of a circular orbit to increase (eq. (B13)), 

As the orbit altitude increases, the half -angle of the cone of tolerance 77 will also in- 
crease (eq. (2)) as the cone of tolerance moves along the ecliptic. In addition, as the or- 
bit altitude increases, the precession rate due to oblateness decreases (eq. (1)). 

Considering all these factors, therefore, it is evident that some relative motion 
exists between the orbit perpendicular and the sunline. This is illustrated in figure 4 
where the orbit perpendicular is shown describing a path within the circle of tolerance. 
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The radius of the circle, which is rj^, is shown increasing with altitude and therefore 
with time. 

An optimum mission (one which maximizes time in continuous sunlight) will t 3 ?pically 
start with the orbit perpendicular on the circle of tolerance at position 1 in figure 4. The 
initial inclination and altitude of the orbit are chosen such that the precession of the orbit 
plane, at first, will be faster than the movement of the Sun along the ecliptic. This 
causes the orbit perpendicular to move across the expanding circle along path a to posi- 
tion 2 where it just approaches the shadow line. At this position, the orbit altitude will 
have increased to the point where the decreased precession of the orbit plane is now 
nearly equal to the movement of the Sun along the ecliptic. Any further increase in alti- 
tude causes the orbit precession to become less than the angular rate of the Sun. This, 
in turn, causes the orbit perpendicular to return across the still expanding circle of 
tolerance until shadow is encountered at position 3. 

For the range of thrust-weight ratios investigated, it was found that the maximum 
time in sunlight is attained when the orbit perpendicular initially starts as far behind the 
sunline as possible which means starting on the edge of the circle of tolerance (such as 
position 1 in fig. 4). Thus, is determined from equation (3) at the initial conditions. 

cos rj„ ft - N sin 6ft cos Ift 

cos ^0 = ^ 0 0 (4) 

N cos 6 q sin Iq 

The angle j^q, therefore, is a function of initial altitude, inclination, and mission start- 
ing date. For a given initial altitude, mission starting date, and thrust -weight ratio, the 
initial inclination must be optimized so that the orbit perpendicular will just approach the 
edge of the circle of tolerance at position 2 in figure 4. An inclination smaller than the 
optimum will cause the orbit perpendicular to follow path b which does not take full ad- 
vantage of the expanding circle of tolerance, thereby reducing the time spent in continu- 
ous sunlight. An inclination greater than the optimum will cause the orbit perpendicular 
to follow path c resulting in a portion of the orbit entering shadow at position 2 on path c. 
Consequently, this will result in a much shorter time spent in continuous sunlight than 
either paths a or b. 

Thrust Reversal 

As discussed in the preceding section, as path a in figure 4 proceeds past position 2, 
the increasing altitude causes the orbit precession to become increasingly less than the 
angular rate of the Sun along the ecliptic. This causes the orbit perpendicular to return 
across the expanding circle of tolerance until, finally, shadow is encountered at posi- 
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tion 3. If, however, the orbit altitude is decreased during the latter portion of the mis- 
sion starting at some point after position 2, the orbit precession can again be increased, 
thereby extending the time spent in continuous sunlight. This can be accomplished by 
simply reversing the thrust vector over the latter part of the mission, thereby reducing 
the altitude. Such a maneuver may be very useful for a development flight such as SERT 
where it is desirable to stay in continuous sunlight as long as possible. 

The effect of thrust reversal is illustrated in figure 5. Without thrust reversal, the 
orbit perpendicular follows path a until shadow is encountered at position 3. With thrust 
reversal, the orbit perpendicular follows path a until the thrust is reversed and thereafter 
follows path b until it finally leaves the shrinking circle of tolerance at position 4. After 
thrust reversal, the altitude, and, therefore, the radius of the circle of tolerance, de- 
creases until it is approximately equal to the initial altitude at position 4. 


Methods of Computation 

Precise solutions for the type of mission described in this report could be obtained 
by numerically integrating each orbit of the trajectory by using a computer code such as 
that described in reference 5. However, this is not a feasible approach for a parametric 
study requiring many solutions because of the computer time required. For example, the 
time for the N-Body code of reference 5 to integrate 1000 circular orbits is nearly an 
hour, and a typical mission consists of at least several thousand orbits. 

Consequently, the approximate procedure developed in appendix B is very useful for 
obtaining a large number of solutions with a minimum of computation time. The equations 
derived in appendix B were programed for solution on a high-speed digital computer. 
Computations were made basically at intervals of 10 orbits; that is, the change over one 
orbit of each of the variables of interest was computed at the beginning of an interval and 
simply multiplied by 10 to get the total change over the interval. Each variable was then 
updated. After investigating several computation intervals, the 10-orbit interval was 
chosen as a good compromise between computational accuracy and computer time re- 
quired. 

The computation of each solution consisted of choosing a specific starting date and 
initial inclination. The right ascension and declination of the Sun were obtained for any 
date and time by using the mean orbital elements of the solar ellipse with respect to the 
Earth from reference 5. The initial angle by which the orbit perpendicular lagged the 
sunline was determined from equation (4) . 

For a specific starting date, each initial inclination was checked to see how long the 
mission could proceed before the orbit entered shadow. This involved an orbit by orbit 
comparison of rj and t]^ as computed from equations (2) and (3). The inclination was 
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then iterated until the optimum value (corresponding to the maximum time in continuous 
sunlight) was found to within at least 0.01°. This process was then repeated for different 
mission starting dates until the optimum starting date was found. 

For the cases where the thrust was reversed during the latter part of the mission, 
the optimum starting inclination was first determined by computing the mission with no 
thrust reversal as discussed previously. This same inclination then, proved to be the 
optimum when the thrust was reversed, since the optimum time of thrust reversal always 
occurred after the first close approach to the shadow line (position 2 in fig, 4) . The time 
of thrust reversal was iterated until the optimum value was found. This same process 
was then repeated for different mission starting dates until the optimum starting date was 
obtained. 


RESULTS AND DISCUSSION 
Cases with no Thrust Reversal 

Figure 6 shows the maximum durations spent in continuous sunlight, the correspond- 
ing final altitudes attained, the corresponding optimum mission starting dates, and initial 
inclinations required as functions of initial thrust-weight ratio. The initial altitudes are 
300, 500, and 800 nautical miles (555.6, 926.0, and 1481.6 km). The results in figure 6 
are presented for the case where the Sun is on the same side of the orbit plane as the or- 
bit perpendicular (i. e. , N = +1) and are characterized by mission starting dates in late 
summer and autumn. The solid curves are the results for thrusting tangentially in the 
plane perpendicular to the sunline. For comparison, the dashed curves display the re- 
sults for thrusting tangentially in the orbit plane. 

For the low thrust-weight ratios studied, thrusting in a plane perpendicular to the 
sunline and, therefore, out of the orbit plane has little effect on the mission duration when 
compared with the case of thrusting tangentially in the orbit plane. For example, in fig- 
ure 5, note that the mission duration for an initial thrust -weight ratio of 5.0x10"® and 
initial altitude of 500 nautical miles (926 km) would be 433 days for thrusting tangentially 
in a plane perpendicular to the sunline compared with 428 days for thrusting tangentially 
in the orbit plane. Also note from figure 5 for these two cases, there is little difference 
in either mission starting date or initial inclination. The required mission starting dates 
for these cases would be September 7 and September 9, respectively, and the required 
initial inclinations would be 107. 5° and 107.9°, respectively. The greatest effect of 
thrusting out of the orbit plane is to reduce the final altitude achieved. This is because 
thrusting out of the orbit plane reduces the in-plane tangential component of thrust and, 
therefore, for a given mission duration, the final altitude attained will be less. In fig- 
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ure 6 for the sample cases, the final altitudes attained would be about 3150 and 3500 
nautical miles (5833.8 and 6482.0 km), respectively. 

The magnitude of the thrust component normal to the orbit plane may approach 
90 percent of the total thrust near the end of the mission (77 « 60°). However, the effect 
of this normal component of thrust on nodal precession is at least an order of magnitude 
less than that due to the oblateness of the Earth and its cumulative effect on orbit inclina- 
tion is to reduce the starting inclination only several degrees by the end of the mission. 
For higher thrust-weight ratios, of course, the effects of thrusting out of the orbit plane 
will be increasingly significant. 

Shown in figure 7 is a typical variation of mission duration or time spent in continu- 
ous sunlight as a function of mission starting date. The curves shown are for an initial 
altitude of 500 nautical miles (926 km) and a thrust -weight ratio of 5.0x10" . Other alti- 
tude and thrust -weight ratio combinations would display similar trends, however. It 
should be noted that each point on these curves is fully optimized; that is, each has an op- 
timum initial inclination and 

For any mission starting date, two solutions, N = +1 and N = -1, are possible so 
that the two curves shown in figure 7 are actually continuous and periodic. However, the 
lower portions of the curves representing solutions at unfavorable mission starting dates 
are not shown. The period from about mid-February to mid-August, for example, is un- 
favorable for N = +1 solutions. Favorable mission starting dates occur in the spring for 
N = -1 cases and in the late summer and autumn for N = +1 cases. 

A qualitative understanding of why a particular mission starting date is more favor- 
able than another may be had by considering the mean yearly variation of the declination 
of the Sun shown in figure 8 and the geometry of sketch (c). In sketch (c), are shown 
edgewise views of typical orbit planes for N = +1 and N = -1, each having an inclination 
greater than 90°. Also shown are the typical relative locations of the orbit perpendicu- 
lars and the sunline insofar as these are determined by inclination and declination at the 
beginning of the flight. In reality, the orbit perpendicular will generally not have the 
same longitude as the sunline. The sunline is shown in the equatorial plane with arrows 
indicating that it will move either above or below the equator as time progresses. 




(The plane of the sketch is the plane containing the sunline and the orbit perpendicular.) 

For an N = +1 case, the most favorable starting dates occur in the autumn such 
that, as time progresses, the Sun moves toward negative declinations. This tends to 
keep the Sun in the vicinity of the orbit perpendicular, thus illuminating the entire orbit 
for a longer period of time. This relative alinment is more important early in the mis- 
sion, when the orbit altitude is small, than later in the mission, when the altitude is 
larger and a greater angle can be permitted between the Sun line and the orbit perpendicu- 
lar. For an N = -1 case, the most favorable starting dates occur in the spring such 
that, as time progresses, the Sun moves toward positive declinations. Note that negative 
declinations are unfavorable for an N = -1 case. 

Cases with Thrust Reversal 

In figure 9 are shown the maximum durations spent in continuous sunlight and the 
corresponding maximum altitudes attained as functions of initial thrust -weight ratio with 
thrust reversal. Also shown in figure 9 are the optimum mission starting dates, initial in- 
clinations, and thrust reversal times. The initial altitudes are 300, 500, and 800 nautical 
miles (555.6, 926.0, and 1481.6 km). The results are presented for N = +1 and for 
thrusting tangentially in a plane perpendicular to the sunline (solid curves) and thrusting 
tangentially in the orbit plane (dashed curves) . Here again, as for the cases with no thrust 
reversal, thrusting out of the orbit plane has little effect on either mission starting date 
or initial inclination. The effects on thrust reversal time and maximum altitude attained 
are somewhat larger. The altitudes presented are maximum and not final, since, with 
thrust reversal, the final altitudes are generally about the same as the initial altitudes. 

By comparing the results presented in figure 9 with those in figure 6, it can be seen 
that the mission durations that can be attained are considerably higher when thrust rever- 
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sal is utilized than for the cases having no thrust reversal. On the other hand, the maxi- 
mum altitudes attained with thrust reversal are considerably less than the final altitudes 
attained with no thrust reversal. For example, at an initial thrust -weight ratio of 
5.0x10“® and an initial altitude of 500 nautical miles (926 km), thrust reversal would 
3 deld a mission duration of 602 days compared with 433 days with no thrust reversal 
(solid curves). The maximum altitude attained would be 2200 nautical miles (4074. 4 km) 
compared with 3150 nautical miles (5833. 8 km) with no thrust reversal. 

Also, by comparing the results presented in figures 9 and 6 , it can be seen that the 

effects of thrust reversal on mission starting date and initial inclination are small. For 

-fi 

example, at an initial thrust-weight ratio of 5. 0x10“ and an initial altitude of 500 nauti- 
cal miles (926 km), the required mission starting date with thrust reversal would be 
August 31 compared with September 7 with no thrust reversal (solid curves). The re- 
quired initial inclination would be 107.2° compared with 107. 5°. It may be noted from 
figure 9 that the time of thrust reversal occurs at about 50 to 55 percent of the total mis- 
sion duration. 

A typical variation of mission duration with mission starting date is shown in fig- 
ure 10 for an initial altitude of 500 nautical miles (926 km) and a thrust -weight ratio of 
5.0x10“ . Here again, each point on these curves is fully optimized, that is, each has an 
optimum initial inclination, time of thrust reversal, and 


Thrust Angle History 

Two typical optimized time histories of r) and are shown in figure 11 for no 
thrust reversal and with thrust reversal, respectively. These curves are presented for 
an initial altitude of 500 nautical miles (926 km), a thrust-weight ratio of 5.0x10“ and 
N = +1. The angle 77 ^ is a function only of altitude (as described by eq. (2)) and repre- 
sents the radius of the circle of tolerance . The angle 77 is the angle between the orbit 
perpendicular and the sunline and is given by equation (3) . The angle 77 may also be 
thought of as the maximum angle in yaw between the thrust vector and the orbit plane dur- 
ing any orbit. The angle between the thrust vector and the orbit plane varies sinusoidally 
over one orbit with the maximum value (equal to 77 ) occurring at the two intersection 
points of the orbit track and the plane perpendicular to the sunline. These intersection 
points are illustrated in figure 2 as points 1 and 3. Without thrust reversal (fig. 11(a)), 
the orbit perpendicular just approaches the edge of the circle of tolerance at about 
225 days and finally leaves the circle of tolerance at about 433 days. With thrust reversal 
(fig. 11 (b)), the orbit perpendicular just approaches the edge of the circle of tolerance at 
about 105 days as the altitude is increasing and again at about 510 days as the altitude is 
decreasing and finally leaves the circle of tolerance at about 602 days. 
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Validity of Approximations 

The N-Body trajectory code described in reference 5 was used to investigate the rela- 
tive effects of thrust, perturbations due to the Sun and the Moon, and the oblateness of 
the Earth on the orbit inclination, node, and eccentricity by integrating a typical mission. 
The secular change of either orbit inclination or nodes due to thrust was an order of mag- 
nitude greater than that due to the Sun and Moon combined. No secular change in orbit 
eccentricity due to either thrust or the Sun and Moon was observed. The effects of the 
Sun and Moon were therefore considered negligible and were validly neglected in the ap- 
proximate solution of appendix B. 

The secular rate of change of the orbit node due to oblateness was approximated by 
considering only the terms through the second harmonic in the potential function (eq, (1)) . 
With this approximation, the secular change of nodes due to oblateness was more than two 
orders of magnitude greater than that due to thrust while the secular change in inclination 
and eccentricity was zero. The change of nodes due to higher order oblateness terms 
(third and fourth harmonics) was determined to be only about 0. 1 percent of that due to the 
total oblateness force. The secular change of either inclination or eccentricity due to 
higher order oblateness terms is very negligible. An insignificant error results, there- 
fore, by neglecting higher order oblateness terms in the approximate solution. 

In order to study further the validity of the approximate solution of appendix B, the 
N-Body code was used to integrate numerically a complete mission and results compared 
with those of the approximate solution. The mission chosen was the overall optimum mis- 
sion with no thrmst reversal for an initial altitude of 500 nautical miles (926 km) and an 

— 

initial thrust-weight ratio of 5.0x10 . The mission starting date was September 6 and 

the initial inclination was 107.61°. Orbit inclination and altitude and the angle by which 
the orbit perpendicular lags the Sun in longitude are shown as functions of the number of 
orbits traversed in figures 12(a), (b), and (c), respectively. It can be seen that the 
changes in these variables, as computed from the approximate solution, are very close to 
the numerically integrated solution obtained with the N-Body code . 

The approximate solution yielded a maximum time in continuous sunlight of 433 days 
and the thrust angle history for this optimum case is shown in figure 11(a). The N-Body 
solution entered shadow at 440 days but also encountered some shadow earlier in the mis- 
sion as illustrated by the thrust angle history shown in figure 11(c). The reason is that 
the N-Body case started with an inclination which was optimum only for the approximate 
solution. For the N-Body case to remain in continuous sunlight, it would be necessary to 
start with an inclination slightly less than that found to be optimum from the approximate 
solution. It is estimated that this required small change in inclination would be about 
0. 1° and would yield a maximum time in continuous sunlight very nearly equal to the 
433 days found with the approximate method. 
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The computer time required for the approximate solution during which 10 complete 
trajectories were flown while iterating on the optimum initial inclination was about 
0. 5 minute. The time required for the N-Body to fly just one trajectory was about 
3 hours. 


Error Tolerances 

The data presented in this report assume that the initial orbit conditions such as the 
initial inclination, altitude, and launch time of day which determines are exact. Fur- 
thermore, it is assumed that the electric rockets on the spacecraft produce full thrust at 
the start of the mission and continue to do so throughout the mission. The values of mis- 
sion duration and altitude attained, thea, are essentially absolute maxima. 

Realistic satellite launchings will always result in some small deviation from the de- 
sired altitude, inclination, and eccentricity. Also, a reasonable launch window is desired 
so that a range of xj/Q (function of launch time of day) must be accounted for as well as er- 
rors in establishing the desired xf/Q. If these deviations in the initial orbit parameters 
are not accounted for, the satellite may begin the mission in shadow or may enter shadow 
much sooner than expected. Dispersions in the initial orbit parameters may be accounted 
for by biasing the initial inclination with only a small penalty in mission duration. 

Biasing the initial inclination was briefly investigated for a typical mission with 
thrust reversal. Typical launch injection errors of ±15 nautical miles (±2.778 km) in or- 
bit altitude and ±0. 5° in orbit inclination were assumed. Errors in orbit eccentricity 
were not considered but would also have to be accounted for in any real launch. 

The effects of biasing the initial inclination may be illustrated from figure 13 which 
shows mission duration as a function of xpQ for lines of constant initial inclination. These 
results are presented for an initial altitude of 485 nautical miles (altitude error of 
-15 n mi (-2.778 km)). This starting condition is more critical than an altitude error of 
±15 nautical miles (±2.778 km), because it results in a smaller (eq. (2)) and, hence, 
a smaller circle of tolerance at the beginning of the mission. The data of figure 13 all 
have optimum thrust reversal times. 

From the initial conditions presented in figure 13, there is only one initial inclination 
and xJ/q which yield the maximum mission duration of 602 days. The required inclination 
at this point is 106.91° and the required xf/Q is 12. 52°. This case takes full advantage of 
the circle of tolerance by starting on the edge of the circle and again just touching the edge 
of the circle at the first close approach as illustrated by positions 1 and 2, respectively, 
in figure 5. Any other combination of initial inclination and i//q within the limit lines 
shown in figure 13 would result in a mission duration (in sunlight) less than the maximum 
of 602 days. For any mission starting date, two conditions limit the possible combina- 
tions of initial inclination and xf>Q; (1) the necessity for starting the mission in sunlight 
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and (2) the requirement to remain in continuous sunlight as long as possible. The limit 
line to the right of figure 13 defines those missions which just start in sunlight on the 
edge of the circle of tolerance but do not touch the circle of tolerance at the first close 
approach. The limit line to the upper left of figure 13 defines those missions which start 
inside the circle of tolerance but just touch the circle of tolerance at the first close 
approach. The intersection of these two lines defines the case having the maximum mis- 
sion duration. This case would just start in sunlight on the edge of the circle of tolerance 
and would also again just touch the circle of tolerance at the first close approach. 

In order to bias the initial inclination properly to account for orbit injection errors, 
it is necessary to choose an inclination such that no ejqjected dispersion in inclination 
from the biased nominal will exceed the limit lines of figure 13 and yet allow a reasonable 
launch window. Orbit altitude errors can be accounted for by choosing the lowest altitude 
expected, which is most critical, and then biasing the inclination at that altitude. This 
was done for figure 13 which is presented for an altitude of 485 nautical miles (898.22 
km), corresponding to the assumed altitude error of -15 nautical miles (-2.778 km). The 
dashed lines shown in figure 13 define the boundaries of initial inclination and re- 
sulting from assuming a launch window of approximately 25 minutes and an inclination er- 
ror of ±0. 5°. Nominally, the launch window would open at point A = 13. 7°) and close 
at point B (if/^ = 7. 5°). The difference in xJ/q of 6.2° is equivalent to 24.8 minutes of 
launch window since the Earth revolves 0.25° per minute. The biased nominal initial in- 
clination is 105. 82° along a line from point A to point B. For this initial inclination, the 
time in continuous sunlight would range from 583 days at point A to 590 days at point B. 

As can be seen, the nominal inclination was biased to 105. 82° so that for an inclination 
error of +0. 5° resulting in an initial inclination of 106. 32°, the starting conditions would 
not exceed the limit lines of figure 13. The critical points are indicated by point C 
(591 days) at the opening of the launch window and point D (598 days) at the closing of the 
window. An inclination error of -0. 5° (initial inclination of 105. 32°) is not critical from 
the standpoint of exceeding the limit lines but will result in considerably less mission 
duration as indicated by points E (498 days) and F (572 days) . 

A launch injection error which has not been discussed as yet but which can affect the 
length of the launch window is the error in establishing the orbit node at precisely the 
proper inertial longitude. This error directly affects but could be accounted for by 
simply moving inside the boundaries on (13.7° and 7. 5°) shown in figure 13, thereby 
reducing the launch window. 

Although the results of figure 13 are presented for an initial altitude of 485 nautical 
miles (898.22 km) (corresponding to an altitude error of -15 n mi (-2.778 km) from a 
nominal of 500 n mi (926 km)), the initial altitude may range from 485 to 515 nautical 
miles (898.22 to 953.78 km) (±15 n mi (±2.778 km)) depending on the accuracy of the 
launch. For altitudes greater than 485 nautical miles (898.22 km), the results would be 


17 





similar to those shown in figure 13 except that the scale on initial inclination would be 
shifted progressively higher for higher altitudes. In other words, the dashed boundary 
of figure 13 would shift progressively lower in the figure with respect to the limit lines. 
This is illustrated in sketches (d) and (e) for initial altitudes of 500 and 515 nautical miles 
(928.00 and 953.78 km), respectively. As seen in sketches (d) and (e), the overall opti- 
mum point at 500 nautical miles (926.78 km) would have an initial inclination of 107.23° 
and at 515 nautical miles (953.78 km) an inclination of 107. 54° compared with 106.91° at 
485 nautical miles (898.22 km). The maximum mission duration for these three altitudes 
is nearly a constant. At the same time, the initial inclination limits for the launch would 
always range from 105. 32° to 106. 32° (lower and upper lines of dashed boundary) deter- 
mined from biasing the initial inclination to 105.82° at 485 nautical miles (898.22 km) 

(fig. 13). The dashed boundary, therefore, will shift progressively lower with respect to 
the overall optimum point as indicated in sketches (d) and (e). 

It can be seen then, that if the initial altitude were 500 nautical miles (926 km) 

(sketch (d)), the mission durations would range from about 594 days at point D down to 
about 462 days at point E. If the initial altitude were 515 nautical miles (953.78 km) 
(sketch (e)), the mission durations would range from about 592 days at point D down to 
about 425 days at point E. These may be compared with 598 days at point D and 498 days 
at point E for starting at an initial altitude of 485 nautical miles (898.22 km). 
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CONCLUDING REMARKS 


It has been shown that a low-thrust, Earth satellite mission planned around the con- 
cept of remaining in continuous sunlight and utilizing tangential thrusting in a plane per- 
pendicular to the sunline may be analyzed by using the approximate numerical method 
presented in this report. This thrust-vector orientation would allow use of a spacecraft 
with fixed solar panels and thrustors, and attitude control could be implemented by using 
horizon scanners and Sun seekers. Although the analysis and results presented in this 
report were obtained for a satellite having a particular thrust vector orientation with re- 
spect to the orbit plane and the Sun, the approach to the problem and method of solution 
would be generally applicable to other low-thrust, circular Earth orbit missions having 
out-of-plane thrusting. The results obtained by using this approximate method exhibit a 
high degree of accuracy when compared with numerically integrated solutions. 

The time that an Earth satellite remains in continuous sunlight can be maximized by a 
choice of proper initial orbit parameters. These initial conditions were determined for a 
range of thrust-weight ratios and initial altitudes for thrusting tangentially in a plane per- 
pendicular to the sunline. For comparison purposes, tangential thrusting in the orbit 
plane was also investigated. For the low thrust-weight ratios studied, the thrust compo- 
nent normal to the orbit plane had little effect on the mission duration. The greatest ef- 
fect of the normal component of thrust was to reduce the maximum altitude achieved. 

The effect of reversing the thrust during the latter part of the mission was also in- 
vestigated. Optimum thrust reversal times occurred at 50 to 55 percent of the total mis- 
sion duration. It was shown that thrust reversal can significantly increase the time spent 
in continuous sunlight. 

All the orbits studied have inclinations greater than 90° so that all orbits are retro- 
grade. Optimum initial inclinations increase with increasing thrust-weight ratio and in- 
creasing initial altitude. The combination of initial altitude, weight of the satellite and 
the initial inclination required is important in determining the orbital payload capability 
required of the launch vehicle. 

Mission duration increases with increasing initial altitude and decreasing initial 
thrust -weight ratio. The maximum altitude achieved, on the other hand, increases with 
both increasing initial altitude and initial thrust -weight ratio. 

Errors in establishing the initial orbit conditions, as well as starting the mission at 
a time other than the optimum, can cause some of the orbits during the mission to en- 
counter shadow. It was shown that shadowing may be eliminated by biasing the initial in- 
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clination to account for these deviations. This would reduce the time spent in continuous 
sunlight below the maximum mission durations presented in this report. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, March 17, 1967, 
124-06-01-02-22. 
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APPENDIX A 


SYMBOLS 


2 2 

A total thrust acceleration, ft/sec (m/sec ) 

B constant (see appendix B) 

C circumferential thrust acceleration, ft/sec^ (m/sec^) 

D constant (see appendix C) 

E constant (see appendix C) 

2 

g acceleration of gravity on equator at surface of Earth, 32.174 ft/sec 

(9.8066352 m/sec^) 

h altitude of satellite above equatorial radius, ft (m) 

I orbit inclination, deg 

J coefficient of second harmonic in gravitational potential function, 

1.62345x10'^ 

K^,K 2 ,Kg,K^ constants (see appendix B) 
constants (see appendix B) 

orbit index equal to +1 if sunUne is on same side of orbit plane as orbit 
perpendicular or equal to -1 if sunline is on opposite side 

Earth equatorial radius, 20925738.0 ft (6378. 165 km) 

radius of satellite from center of Earth, ft(m) 


K* ,K** 
N 


% 


t 

u 


W 

X,Y,Z 
X’,Y’, Z’ 


time, sec 

instantaneous angular position of satellite in orbit plane measured from 
ascending node positively in direction of motion, deg 

thrust acceleration normal to orbit plane measured positively in the di- 

2 2 

rection of the orbital angular velocity vector, ft/sec (m/sec ) 

Earth centered inertial rectangular coordinates with X pointing in direc- 
tion of vernal equinox and with X and Y lying in equatorial plane and 
Z in right-hand relation 

Earth centered instantaneous rectangular coordinates referenced to po- 
sition of Sun with Z’ pointing in direction of Sun and X’ perpendicu- 
lar to Z’ and lying in equatorial plane and Y’ in right-hand relation; 
primes denote coordinate system used when Sun is in on side of orbit 
plane opposite to orbit perpendicular 
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II 



X.r^Y’',Z" 


a 

r 

5 

V 

% 

9 


M 

O) 

cQd 

Subscripts : 
av 
max 
min 
0 
1 

Superscript: 

♦ 


similar to X’, Y', Z' except used for case where Sun is on same side of 
orbit plane as orbit perpendicular 

right ascension of Sun measured positively counterclockwise from vernal 
equinox to projection of sunline in equatorial plane, deg 

angle by which 6 leads (or lags) u, 0-Nu, deg 

declination of Sun measured positively above equator, negative below, deg 

angle between orbit plane and plane perpendicular to sunline or between 
orbital angular velocity vector (or its negative whichever yields smaller 
angle), and sunline, deg 

maximum t] allowable for orbit to remain in continuous sunlight, deg 

angular position in plane perpendicular to sunline (corresponds to u in or- 
bit plane) measured positively from X' toward Y’ counterclockwise 
about Z’ (sunline opposite to orbit perpendicular) or measured positively 
from X" toward Y" counterclockwise about Z” (sunline on same side 
as orbit perpendicular) ; angle 9 is obtained by projecting satellite posi- 
tion in the orbit into plane perpendicular to sunline along great circle 
which forms spherical isosceles triangle with orbit plane and plane per- 
pendicular to sunline, deg 

gravitational constant of Earth, 1.4076539x10^® ftVsec^ (3.9860319x10^^ 
m Vsec^) 

angle by which orbit perpendicular lags Sun in longitude, a - eil>, deg 

angular velocity of satellite, rad/sec 

right ascension of orbit perpendicular, measured positively counterclock- 
wise from vernal equinox to projection of orbit perpendicular in equatorial 
plane, deg 


average 

maximum 

minimum 

initial conditions 

intermediate values for matrix 

indicates constant angle for one orbit (see appendix C) 
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APPENDIX B 


APPROXIMATE NUMERICAL SOLUTION FOR LOW THRUST ORBIT RAISING FROM 
CIRCULAR ORBIT WHILE REMAINING IN CONTINUOUS SUNLIGHT 
AND THRUSTING TANGENTIALLY IN PLANE 
PERPENDICULAR TO SUNLINE 

This appendix derives an approximate numerical solution for the normal component of 
thrust acceleration W and then determines its effects on orbit inclination and ascending 
node. These effects can then be combined with those due to the oblateness of the Earth. 

Consider the geometry of a satellite in a retrograde (inclination greater than 90°) 
circular orbit about the Earth. The orbital angular velocity vector will be defined as the 
orbit perpendicular. For generality, the Sun may be either on the same side of the orbit 
plane as the orbit perpendicular (Z" in sketch (f)) or on the opposite side (Z’ in 
sketch (f)). 
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Since the thrust vector is constrained to be always tangential and in a plane perpen- 
dicular to the sunline, the acceleration components perpendicular to the simline and along 
the sunline are 


A^, = Ay, = -NA sin 9 
Ay» = Ay,, = NA cos 9 > 

A^i — ^2” ~ ^ 


(Bl) 


where 9 is measured differently between the X' and X" systems (see appendix A) . 

These acceleration components can now be expressed in the X, Y, Z system by ro- 
tating the X’,Y',Z* system into the X, Y, Z system. 

For the case where the sunline is on the side opposite to the orbit perpendicular, the 
X’, Y’, Z’ system may be rotated counterclockwise about X' by an amount equal to 
(90° - 6). 



0 0 


X’ 

sin 6 cos 6 


Y* 

cos 6 sin 6 


Z’ 


(B2) 


For the case where the sunline is on the same side as the orbit perpendicular, the 
X", Y”, Z" system can be rotated counterclockwise about X" by an amount equal to 
[90° + (-6)] = (90° - 6). This rotation (N = +1) is identical to the preceding case (N = -1), 
and equation (B2) then describes the first rotation for both N = +1 and N = -1 . 

The second rotation is needed now to aline the Xj, Yj, Z!^ system with the X, Y, Z 
system . 

For the (N = -1) case, the Xp Y^, Z'^ system can be rotated clockwise about Zj by 
an amount equal to -(a - 90°) = (90° - a) or, for the (N = +1) case, the Xj, Yj, Z'^ sys- 
tems can be rotated clockwise about Zj by an amount equal to -(a - 270°) = (270° - a). 

Here again, the rotation for either case is equivalent and can be written as 


X 


sin or cos or 0 



Y 

= 

-cos or sin a 0 


Yi 

Z 


0 0 1 



- _ 


- 




(B3) 
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Now, combining the two rotations by matrix multiplication yields 


X 


sin a cos a 0 


10 0 


X’ 

Y 

= 

-cos a sin a 0 


0 sin 6 cos 6 


Y’ 

Z 


0 0 1 


0 -cos 6 sin 6 


Z» 

X 


sin a cos a sin 6 

cos a cos 6 


X’ 

Y 

= 

-cos O' sin a sin 5 

sin a cos 6 


Y' 

Z 


0 -cos 6 

sin 5 


Z» 


As is explained later, the rotation matrix of equation (B4) can be considered 
Therefore, twice differentiating equations (B4) yields 


= Aj^, sin a + Ay> cos a sin 6 + A^i cos at cos 6 
Ay = cos a + Ay, sin a sin 6 + Ag,, sin o; cos 6 

'Z’ 

Now, substituting from equations (Bl) yields 


A^ = 0 - Ay, cos 6 + A.7, sin 6 


Aj^ = -NA sin 9 sin a + NA cos 9 cos a sin 6 


Ay = NA sin 9 cos a + NA cos 9 sin a sin 6 




A^ = -NA cos 9 cos 6 


Since 9 is an unknown angle for this problem, it is now necessary to define 
of other known angles. From spherical trigonometry (see appendix C): 

sin 9 = DN sin u + E cos u 

cos 0 = D cos u - EN sin u 


Therefore, by substituting equations (B8) into equations (B5) 


(B4) 

constant. 


(B5) 


9 in terms 


(B6) 
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= -Na[(DN sin u + E cos u)sin O' - (D cos u - EN sin u)cos a sin 
Ay = NA[(DN sin u + E cos u)cos ot + (D cos u - EN sin u)sin a sin 6^ 
A^ = -Na[(D cos u - EN sin u)cos ^ 


or 


A^= -Na[(DN sin a. + EN cos a. sin 6 )sin u - (D cos a sin 6 - E sin d)cos 
Ay = NA[(DN cos Q!-EN sin a. sin 6)sinu + (D sin a sin 6 +E cos 0!)cos 

EN cos 6 sin u + D cos 6 cos ul 



y (B7) 


Now, from reference 6, the perturbative acceleration due to the thrust component normal 
to the orbit plane is 


W = A^ cos cfis sin I + Ay sin sin I + A^ cos I (B8) 

and is measured positively in the direction of the angular velocity vector. Note that the 
angle as used in this development is 270° greater than the angle cQ> from refer- 
ence 6. Substituting equations (B7) into equation (B8) yields: 

W = NA {[-(DN sin a + EN cos a sin 6)sin u + (D cos a. sin 5 - E sin 0 !)cos v^cos cQ, sin I 

+ [(DN COS Of - EN sin a sin 6)sin u + (D sin a sin 6 
+ E cos Q!)cos ujsin cQo sin I - [^EN cos 6 sin u +D cos 6 cos i^cos l| 
Now multiplying and combining terms yields 
W = NA|dN cos a sin Sh sin I sin u - DN sin a cos sin I sin u 

- EN sin Of sin 6 sin cfia sin I sin u - EN cos a sin 6 cos sin I sin u 
+ EN cos 6 cos I sin u + D cos a sin 6 cos <£h> sin I cos u 
+ D sin a sin 6 sin sin I cos u - E sin a cos cfia sin I cos u 
+ E cos a. sin <Slo sin I cos u - D cos 6 cos I cos 
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or 


W = Na|dN sin I sin u(cos a sin efio - sin a cos ciL) - EN sin 6 sin I sin u(sin a sin 

+ cos a cos oJii) + EN cos 6 cos I sin u - D cos 6 cos I cos u 
+ D sin 6 sin I cos u(cos a cos <£U, + sin a sin cQo) 

- E sin I cos u(sin a cos cfio - cos a sin 

Now, 


\p = a - t£h 


Therefore 

W = Na[(-DN sin I sin i// - EN sin 6 sin I cos + EN cos 6 cos I)sin u 

+ (D sin 6 sin I cos - E sin I sin xp - D cos 6 cos I)cos 


Also, from reference 6, for a circular orbit 


cQd = 



1^0 


(B9) 


(BIO) 


In order to arrive at a closed form solution for I and efis, it is necessary to assume 
that over the period of one orbit, none of the orbit parameters change, except u, of 
course; that is, for one orbit, I, <Sh>, a, 6, \f/, and r are constant with u being the in- 
dependent variable of integration; therefore, let 


B = 



Kj^ = D (see appendix C for D and E) 
K2=E 
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Kg = sin I sin x}/ 


= sin 6 sin I cos - cos 5 cos I 

Then, W can be written as 

W = -AjfKgK^ + KjKg)sin u + N(K2Kg - KjK4)cos uj (Bll) 

Substituting equation (Bll) into equations (B9) and (BIO), yields 

i = -Ab|(K2K4 + KjKg)sin u cos u - N(KjK4 - K2Kg)cos^ 

^ + KjKg)sin^ u - N(K^K^ - K2Kg)sin u cos 


Now, let 


Then, 


K* = -AB(K2K4 + KjKg) 
K* * = NAB(KjK 4 - K2Kg) 


I = K* sin u cos u + K* * cos 


u 


• 1 2 

efio = (K* sin u + K* * sin u cos u) 

sin I 

Since u is the independent variable, integration must be made with respect to u. Let 

du = w dt 


where co is the angular velocity of the vehicle and is assumed to be constant over one 
orbit . Therefore , 

1 2 
dl = — (K* sin u cos u + K* * cos u)du 

u> 


dcf2o = — (K* sin^ u + K* * sin u cos u)du 


w sin I 
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Integrating yields 



The periodic terms over one orbit, can be neglected and the secular changes in I and cil 
due to W are just 


AI = — K* * Au 
2co 


A cilo — 


2u> sin I 


K* Au 


(B12) 


From reference 2, the change of altitude as a function of time for a circular orbit is 
given by 


dt Rg g ^ 


+ h) 


3/2 


Now, since 


du = CO dt 


— v£ — dt 
(Rj, + h)3/2 


then 


dh = — (Rg + h)^ du 

Vm Rg Vg 


(B13) 


The circumferential thrust acceleration C is a function of u and varies periodically 
over one orbit. Therefore, over one orbit, C can be taken as a constant, that is, 

C + C 

Q _ max min 
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where 


and 

Integrating equation (B13) 3 delds 

P p 40 

(Rg + hg) - (Rg + h) "^ = ^ An (B14) 

Vm Rg Vg 

From equation (B14), the change in altitude from hg to h can thus be found with a 
change in u, Au. 

Also from reference 2, the nodal precession of the orbit due to oblateness can be ex- 
pressed as 


cd, « - Vm JRg^(Rj, + cos I 


therefore, 




o)R 


3/2 


E 


VRg + hj 


Substituting for co and for du from equation (B13) gives 

dSL = - 7 ^ cos I(Rg + h)"® dh 


20 


av 


Integrating 


Acft. = - Vitg ^ cos ijcRg + Iiq)"'* - (Rg + h) ^ 


(B15) 


Therefore, from equation (B15) the nodal precession due to oblateness can be found after 
the change in altitude has been found. This can be combined with equation (B12) to give 
the total change in nodes. 
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APPENDIX C 


EQUATIONS DEFINING ANGULAR POSITION IN PLANE PERPENDICULAR 
TO SUNLINE IN TERMS OF KNOWN ANGLES 

Consider the spherical triangle formed by the orbit plane, the plane perpendicular to 
the sunline and the equatorial plane in sketch (g) . Let u* , 9* , and y* be angles which 
have a long term variation with time but may be considered constants for one or even 
several orbits. The angle u* is measured positively from the ascending node (u = 0) to 
the point of intersection of the orbit plane with the plane perpendicular to the sunline; 9* 
is measured negatively from the equator to the intersection of the planes; and y* is 
measured with the same sense as 9* and is the amount that 9 leads (or lags) u (see 



appendix A for definitions of u and 9). When u = 0, 9 =y* . Now, from sketch (g). 


Nu* =9* - y* 


for 


0°<u* <+180° 
-180° <0* < 0° 


or 
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y* = Q* . NU* 


Also, 


9 = y* + Nu + 360° = (0* - Nu* ) + (Nu + 360°) 


Therefore, expanding yields 

sin 9 = sin(0* - Nu*)cos(Nu + 360°) + cos(0* - Nu*)sin(Nu + 360°) 
sin 9 = (sin 9* cos Nu* - cos 9* sin Nu* ) (cos Nu cos 360° - sin Nu sin 360°) 


+ (cos 9* cos Nu* + sin 9* sin Nu* ) (sin Nu cos 360° + cos Nu sin 360°) 


Now, since 


cos Nu* = cos u* 


and 


sin Nu* = N sin u* 

sin 9 = (sin 9* cos u* - N cos 9* sin u* )cos u + (cos 9* cos u* + N sin 9* sin u* )N sin u 
or 

sin 0 = N sin u(cos u* cos 0* + N sin u* sin 0* ) 

+ cos u(cos u* sin 0* - N sin u* cos 0*) (Cl) 


Similarly, 

cos 0 = cos(0* - Nu*)cos(Nu + 360°) - sin(0* - Nu*)sin(Nu + 360°) 
cos 0 = (cos 0* cos Nu* + sin 0* sin Nu* ) (cos Nu cos 360° - sin Nu sin 360°) 

- (sin 0* cos Nu* - cos 0* sin Nu* ) (sin Nu cos 360° + cos Nu sin 360°) 

Therefore, 
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cos d = cos u(cos u* cos 0* + N sin u* sin 9 * ) 


- N sin u(cos u* sin 0* - N sin u* cos 0* ) (C2) 

From equations (Cl) and (C2), let 

D = cos u* cos 0* + N sin u* sin 9* 

E = cos u* sin 0 * - N sin u* cos 0* 


Therefore, 


sin 0 = DN sin u + E cos u 
cos 0 = D cos u - EN sin u 

Now u* and 0* can be found in terms of other known angles. Referring again to 
sketch (g) and using the law of sines and the law of cosines gives, for N = -1, 

sin u* = sin(90° - 5)sin(;^ + 180^ _ _ cos 5 sin \f/ 
sin T} sin tj 

cos u* = - 5) + cos I cos ri _ sin 6 + cos 77 cos I 

sin I sin 77 sin I sin r) 


and, for N = +1, 


sinu* - + 6 )sin _ 

sin 7] 




sin 77 


cos 


_ c os(90° + 5) + cos 1 cos 77 _ -sin 6 + cos I cos 77 
sin I sin 77 sin I sin 77 


Therefore, in general, 

^ TWT COS 6 sin xL/ 
sin u* = N = — ^ 

sin 7 ] 

> 

cos u* = -N sin 5 + cos 7 ? cos I 

sin 77 sin I ^ 


(C3) 
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Likewise, for N = -1, 


sin 9* = 


_ sin I sin(;^ + 180°) _ sin I sin rj/ 


sin 


sin rj 


cos 9* = ^ + cos(90° - 6)cos rj _ cos I + sin 6 cos r] 

sin(90° - 6) sin 77 cos 6 sin 77 


and, for N = +1, 


sin 


g* _ _ sin(180° - I)sin if/ _ _ sin I sin ^ 
sin 77 sin rj 


cos 0* = 


_ cos(180° - I) + cos (90° - 6)cos rj _ -cos I + sin 5 cos rj 


sin(90° - 6) sin rj 


cos 6 sin rj 


Therefore, in general. 


sin 0* = - N sin I sin x}y 
sin tj 

cos 0* = cos I + sin 5 cos rj 
cos 6 sin 77 





(C4) 


It is seen that u* and 0* are functions only of 6, rj, and I, which are assumed 

constant for one orbit. In order for equations (C3) and (C4) to be evaluated correctly, it 
is important that, for N = +1, and for N = -1, if/ = xf/ + 180°. 
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Figure 1. - Conceptual spacecraft configuration with solar cell panels. 
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Ang^e by which orbit perpendicular lags sun in longitude, ip, deg 

Figure 4. - Typical paths of orbit perpendicular within circle of tolerance. 
Orbit index, +1. (Letters denote p^h of orbit perpendicular; numbers 
denote position of orbit perpendicular. ) 
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Angle by which orbit perpendicular lags sun in longitude, vp, deg 


Figure 5. - Typical optimum paths of orbit perpendicular within circle of tolerance with and 
without thrust reversal. Orbit index, +1, 





Declination of sun, 6, deg Mission duration, days in sunlight 



Mission starting date 


Figure 7. - Maximum mission duration as function of mission starting date. Initial altitude, 
500 nautical miles; initial thrust-weight ratio, 5.0x10"^; tangential thrusting in plane 
perpendicular to sunline. 



Figure 8. - Mean declination of sun as function of mission starting date. 
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Time of thrust reversal, days Mission duration, days in sunlight 







Mission duration, days in sunlight 



Mission starting date 


Figure 10. - Mission duration as function of mission starting date with thrust reversal. Initial 
altitude, 500 nautical miles; initial thrust-weight ratio, 5. 0x10" ; tangential thrusting in 
plane perpendicular to sunline. 
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(cIN'lxxly solution with no thrust reversal. Mission starting date, September 6; initial 
inclination, 107.61° initial angle by which orbit perpendicular lags sun in longitude 
16.04°. 

Figure 11. - Comparison of angle between sunline and orbit plane and maximum allowable 
angle for orbit in continuous sunlight for typical missions. Initial altitude, 500 nautical 
miles; initial thrust-weight ratio, 5.0x10“^; orbit Index, +1. 


Orbit altitude, n mi QpDft rncllnation, deg 



(a) Orbit inclination as function of number of orbits traversed for typical mission in continuous 
sunlight. 



(b) Orbit altitude as function of number of orbits traversed for typical mission in continuous 
sunlight. 

Figure 12. - Comparison of approximate solutions with N-body solution both as function of number 
of orbits traversed. No thrust reversal; Initial altitude, 500 nautical miles; Initial thrust-weight 
ratio, 5.0x10"^; orbit index, +1; mission starting date, September 6; Initial angle by which orbit 
perpendicular lags sun in longitude, 16.04®; initial inclination, 107.61°. 
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(c) Angle by which orbit perpendicular lags sun in longitude as function of number of orbits 
traversed for typical mission in continuous sunlight. 

Figure 12. - Concluded, 



lags sun in longitude, deg 

Figure 13. - Mission duration as function of initial angle by 
which orbit perpendicular lags sun in longitude. Thrust 
reversal; initial altitude, 485 nautical miles; initial thrust- 
weight ratio, 5.0x10"^; mission starting date, August 31; 
orbit index, +1. (Letters denote discrete points at opening 
and closing of launch window. ) 
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aeronautical and space activities of the United States shall he 
conducted so as to contribute to the expansion of human knowU 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof/* 

— ^National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 
importance as a contribution to existing knowledge, 

TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Scientific and technical information generated 

under a NASA contract or grant and considered an important contribution to 
existing knowledge. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 
activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
non-aerospace applications. Publications include Tech Briefs, Technology 
Utilkation Reports and Notes, and Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

WashinstoH/ D.C. 20546 


